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■ Creation of entanglement is considered theoretically and numerically in an ensem- 



(N 



O 



ble of spin chains with dipole-dipole interaction between the spins in each chain. The 
unwanted effect of the long-range dipole interaction is compensated by the optimal 
choice of the parameters of radio-frequency pulses implementing the protocol. The 



>^ • errors caused by (i) the influence of the environment, (ii) non-selective excitations, 

, (iii) influence of different spin chains on each other, (iv) displacements of qubits 

ff^ ' from their perfect locations, and (v) fluctuations of the external magnetic field are 

o ■ 

^ ' estimated analytically and calculated numerically. For the perfectly entangled state 



the z component, M, of the magnetization of the whole system is equal to zero. 
The errors lead to a finite value of M . If the number of qubits in the system is 



' large, M can be detected experimentally. Using the fact that M depends differently 

O^' on the parameters of the system for each kind of error, varying these parameters 



would allow one to experimentally determine the most significant source of errors 



H ' and to optimize correspondingly the quantum computer design in order to decrease 

the errors and M. Using our approach one can benchmark the quantum computer, 
decrease the errors, and prepare the quantum computer for implementation of more 
complex quantum algorithms. 

PACS numbers: 03.67.Lx, 75.10.Jm 



I. INTRODUCTION 



A quantum computer is supposed to be an analog device designed for implementation 
of quantum algorithms. The most important algorithm is the Shor's quantum algorithm 
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for factorization of large integer numbers Before experimental realizatation of this 

and other complex quantum algorithms, this computer must be tested by implementation of 
simple algorithms, such as, for example, creation of entanglement with two quantum states 
(and many qubits). This procedure would allow one to identify the most significant sources 
of errors and to optimize the quantum computer design. 

The following criteria must be met by a physical implementation of a quantum infor- 
mation processor: (i) the qubits must be easy to physically manipulate, (ii) easy to increase 
the number of qubits, (iii) qubits should interact with each other, (iv) the qubits should 
be somewhat isolated from environment, and (v) the qubit states must be detectable. Cur- 
rently, solid state implementations of a quantum computer allow (theoretical) scalability 
[criterion (ii)] and possibility to substantially increase the number of qubits in the quantum 
computer register. In a solid state currently one and two-qubit quantum logic operations are 
implemented and measured experimentally. Nuclear magnetic resonance and electron spin 
resonance techniques have been used for quantum information processing in solids. Etan- 
glement between electron and nuclear spins of the same molecule has been implemented 
and measured in a malonic acid single crystal ^| and in ^^N^lCeo endohedral fullerene ^j. 
Quantum process tomography has been performed on a solid-state qubit represented by a 
nitrogen-vacancy defect in diamond ^6] at room temperature. Initialization (cooling) of nu- 
clear spin has been realized in isotopically labeled malonic acid molecules by a controllable 
transfer of polarization from neighboring nuclear spins 0,01 at room temperature. One and 
two-qubit quantum logic operations have been implemented and measured in superconduct- 



ing quantum computers jSl, US 



A further technological advance can be implementation 



of simple quantum computing algorithms in a scalable solid state system with a large (more 
than two) number of qubits for demonstration of basic principles of quantum computation. 
In this paper, we consider how to implement, probably, the simplest possible algorithm with 
many qubits: creation of entanglement. To implement this algorithm, we choose one of the 
most affordable systems: a spin chain placed in a permanent magnetic field with a gradient 
along the chain. There is a constant dipole-dipole interaction between the qubits. The logic 
operations are implemented by rectangular radio-frequency pulses. This setup allows one to 
achieve all criteria (i)-(v). 

Indeed, criterion (i) is satisfied because in the system there are no switchable interactions 
controlled by nanoscopic metal gates, and quantum logic operations are implemented using 
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global addressing technique based on rectangular radio-frequency pulses. Criterion (ii) is 
satisfied because we consider solid-state QC architectures and because rectangular pulses 
can be used in the QC with many qubits and all parameters of the applied pulses can be 
determined analytically. The number of pulses in our entanglement protocol is equal to the 
number of qubits in the chain. The long-range constant magnetic dipole-dipole interaction 
satisfies criterion (iii). In the majority of spin-based quantum computer architectures this is 
the only interaction when neighboring qubits are separated from each other at the distances 
> 1 nm. Criterion (iv) is satisfied because our system allows implementations using qubits 
with long decoherence times. One possible implementation is based on nuclear spins, for 
example nuclear spins of ^^P ^^Si Q, or Li Q in ^^Si. Another implementation is 
based on endohedral fullerenes, ^^N@C6o and ^^PQCgo [l^, where the fuUeren cage provides 
a good isolation of electron or nuclear spin of the enclosed atom of nitrogen or phosphorus 
from the environment. The electrons of the nitrogen or phosphorus atoms have spin 3/2, and 
the nuclei have spin 1/2 1^, Consequently, our protocols can work only for the nuclear 
spins of the endohedrals. However, a modification of our scheme for the electrons with spin 
3/2 is possible The third possible implementation is based on electron spins in self- 



assembled monolayer systems 
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2l|. Measurement [criterion (v)] can be performed 
by using an ensemble of many identical spin chains shown in Fig. 1. The external 
magnetic field in Fig. 1 is nonuniform in the x direction but uniform in the y direction. 
Such a gradient can be created by current (s) flowing along the y axis in the plain formed by 
the chains. Consequently, we take the angle G between the qubit plane and the permanent 
magnetic field to be B = 7r/2. If the distance between neighboring chains is sufficiently 
large all chains experience the same conditions. Consequently, one can create entanglement 
simultaneously in all chains. If the number of spins in the system is sufficiently large, one 
can detect the macroscopic magnetization Ai of the whole system. If each chain is in the 
entangled state 

-^(|00...00) + e^111...11)), (1) 

then z component Ai^ of the macroscopic magnetization must be equal to zero. All possible 
errors result in deviation of from its perfect value Aiz = 0. We will show in this paper 
that the sign of and the dependence of on different parameters of the model can 
provide us the information about the main source of error. This information can be useful for 
benchmarking the quantum computer and optimization of its architecture and parameters 
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B" (arb. units) 
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FIG. 1: A setup for creation of entanglement with ensemble of spin chains. 



of quantum protocols. 



II. HAMILTONIAN 



We consider the implementation of entanglement only for one spin chain of the ensemble 
of chains illustrated in Fig. 1. Since the chains experience the same conditions the dynamics 
of all chains is similar. One can neglect the magnetic dipole-dipole interaction between the 
spins of the different chains when the condition (a/d)^ <^ 1 is satisfied. (Here a is the 
distance between neighboring qubits of a single chain and d is the distance between different 
chains, see Fig. 1.) More exact condition of relative independence of different chains will be 
given in Section FV (^1 below. 

The applied magnetic field has the following components 

Bn{x, t) = Bl cos(z/„)f: + ipn), -Bl sm{unt + <^n] 

where B^{x) is a permanent magnetic field and 5^, Un, and ipn are, respectively, the 
amplitude, frequency, and phase of the nth circularly polarized radio-frequency rectan- 
gular pulse of a protocol. The permanent magnetic field B^{x) has a constant gradient 
B^{xi) = ujq + 16 B^, where xi is the x-coordinate of the /th qubit. Note that in practice the 
magnetic field gradient can be variable. Our results can be easily reformulated for this case. 
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provided that — ^ B^^. If there is the dipole-dipole interaction between 

the qubits the Hamiltonian is 

H^ = H^ + H,^, + Vn{t), (2) 

where 

H^ = -J2 ^iSt, Vn{t) = [ST exp [-1 {i^nt + <^„,)] + h.c] , (3) 

= (^Sf5|. J=^(3cos=e-l) = -^, (4) 

Here the Hamiltonian is presented in the frequency units; h is the Planck constant; Sf^ = 
Sf ± iSf; Sf, Sf , and are the components of the operator of the fcth nuclear or electron 
spin 1/2; Ui = Uq + I6u!, Uq = 7_B°(xo), Su = ■ySB^; Qn = 7-8^ Rabi frequency of the 

rath pulse; 7 = or 7 = — 7e, where •yiy and 7e are, respectively, the nuclear and electronic 
gyromagnetic ratios; A = 1 nm, A is the dimensionless parameter equal to the distance 
between neighboring qubits measured in nanometers, so that a = AA, (a is the distance 
between neighboring qubits); G = 7r/2 is the angle between direction of the spin chain and 
direction of the permanent magnetic field; /i = — (^e/^B for an electron spin and /x = QnI^n 
for a nuclear spin, /xb and are, respectively, the Bohr and nuclear magnetons, Qe ^ 2 
and gjy are, respectively, the electron and nuclear g-iactors. In Eq. ^ we neglect the x 
and y components of the dipole-dipole interaction because the magnetic dipole field on Ith 
qubit in any stationary state |00 . . . 00), |00 . . . 01), ... is much less than 6B^. So, only the 
z component of the dipole-dipole field significantly affects the quantum dynamics. 



III. CREATION OF ENTANGLEMENT 



First, we discuss the formal steps required to create the entangled state. Let the initial 
state be the ground state \0l-i0l-2 ■ ■ ■ OiOq) as shown in Fig. 1. We split the ground state 
into two states by applying the Hadamard transformation 

^o|Ol-iOl-2 . . . OlOo) = -^(|0l_i0l_2 . . . OiOo) + e'^' \0l-i0l-2 ■ ■ ■ Oilo)). (5) 

Here and below we omit the total phase factor. Then we apply Control-Not gate CNqi to 
the 1st qubit to obtain 

-^i\0L-i0L-2 ■ . ■ OlOo) + e^^"|0i_i0i_2 . . . 02lilo)). (6) 
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Applying Control-Not gates to the remaining L — 2 qubits we obtain the entangled state: 

-^(|0i-i0z.-2 . . . OlOo) + . . . lilo))- (7) 

Below we do not take into consideration the phase factor e**^. One can make to be equal 
to zero by a proper choice of the phases of the pulses js^. 



IV. QUANTUM DYNAMICS 



We decompose the wave function into the basis states \p) of the unperturbed Hamiltonian 



2^-1 



m = E C,{t)e-^^^'\p), 
where \p) = \nL_1nL-2 ■ ■ .rii. . . niUo), rii = 0, 1, and 



L-l 



E, = {p\H' + H,M = - E -iip>Kp) - ^ E _ E^ 7^^- 



(8) 



(9) 



1=0 1=0 k=l+l 

Here sf{p) is the eigenvalue of the operator S'f in the state \p): if the Ith spin of the state 
\p) is in the state | . . . 0; . . .) then sf{p) = 1/2 and if the Zth spin of the state \p) is in the 
state I . . . 1; . . .) then sf{p) = —1/2. 

Let the jth spin of the state \p) be in the state \ ...0j...) and let the state \q) be 
associated with the state \p) by flip of the jth spin, \q) = \ . . .Ij . . .). If the frequency z/„ 
of the electromagnetic field is close to the Larmor frequency of the jth spin, then under 
the conditions I^cjI ^ JM^ and |5u;| ^ Qn the pulse affects mostly this jth spin, and the 
approximate solution is j23| 



COS 



K{q,P)Tn 



, ■^n{q,p) 
+ i- — : r sm 



Cq(tn + T„) 



t— -sm 

K{q,P) 

Cp(tn) = 



^n{q,p) 
A„(g,p)r„ 



>^n{q,p)rn 



-iA„{q,p)T„/2 



JAn{q,p)tn-iiPn+iAniq,p)rn/2 



(10) 



1, Cg{tn) = 0, 

where t„ is the time of the beginning of the nth pulse and 



Aniq.p) = Eq- Ep- A„(g,p) = ,J Al{q,p) + ^1. 

If the detuning A„(g,p) is equal to zero, A„(g,p) = 0, and for flnTn = (7r-pulse), there 
is a complete transition between the states \p) and |g). Here we neglected the transitions 
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associated with flips of spins witli nonresonant transition frequencies for which I ^ j. The 
corrections to the probabihty amphtude associated with flips of these spins are of order of 



r2„/(2|5ci;| |j — /|) [2J,|2^. These corrections are small provided f2„/|5cc;| -C 1. 



A. Control-Not gate for the system with long-range dipole-dipole interaction 

Since the entanglement is implemented by a sequence of Control-Not gates we now derive 
the parameters required to implement these gates taking into consideration the long-range 
interaction. Because of this interaction the action of the Control-Not gate on the jth qubit 
depends not only on the states of (j — l)th and (j + l)th qubits, but also on the states of 
distant Ith qubits, where / 7^ j — 1, j, j + 1. Consequently, it is convenient to formulate the 
Control-Not gate CN not in terms of the states of the corresponding qubits, like CNj_i j, 
but in terms of the eigenstates \p) of the Hamiltonian H'^. For example, the gate CN j{p, q) 
flips the jth qubit for the state \p) and completely suppresses flip of the same jth qubit for 
the state |g). This procedure allows one to compensate the unwanted action of the long- 
range interaction by an optimal choice of the parameters of the pulses. Our approach works 
for the case when there are only two "working states" (excluding error unwanted states) in 
the quantum register. It is useful for creation of entanglement with two states introduced 
in Section IIIII If there are more than two states in the quantum register, more complex 
(shaped) pulses are required to minimize the effect of the long-range interaction. 

Let us write the entangled state in the form (l/-\/2)(|0) + |p')), where |0) is the ground 
state and is the excited state. The pulses of the protocol must suppress the transitions 
from the ground state and to implement the transitions for the excited states, so that the 
state evolves as 

\p') : |000 . . . 001) ^ |000 . . . Oil) ^ |000 ... 0111) ^ ... ^ 

|001 . . . Ill) |011 . . . 11) ^ |111 . . . 11). (11) 

The Control-Not gates implementing these transitions are CNj(p', 0), where j = 1, 2, . . . , L — 
1 is the number of the spin to be flipped in the excited state \p') and also the number of the 
corresponding vr-pulse. (We do not count the initial 7r/2 pulse.) 

We now derive the parameters of the initial pulse implementing the Hadamard transform. 
In order to flip the 0th spin in the ground state the detuning must be equal to zero, Ei — 
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Eq — vh = 0, where vh is the frequency of the pulse. Using Eq. Q we obtain 

i^H = Er-Eo = u;o + ^Y.Jr (12) 

The Rabi frequency VLh of the initial pulse must satisfy the condition VLh ^ and the 
time-duration of the pulse is th = ii/{2VLh)- 

Consider now the Control-Not gate CNj(p',0). If the state \q') is associated with the 
state \p') by a flip of the jth spin, the frequency is 

^, = E,-E, = u:, + -Y.^^. (13) 

For example, for j = 1, one has Sq{p') = —1/2 and sf{p') = 1/2 for / = 2, 3, . . . , L — 1, so 
that the frequency of the flrst vr pulse is 

J 1 

1=3 

the frequency of the second n pulse is 



J 1 
= + 2l^ g (7^' 

and so on. 

In order to suppress the unwanted transitions from the ground states one can use a 
2itK method 0, ■ Here we will modify this method so that it can be applied to the 

system with the long-range interaction. As follows from Eq. |T0|) the transition with nonzero 
detuning is suppressed if the value of the sine is equal to zero, i.e. when 

where = 1, 2, . . . and the state \q") is associated with the ground state by a flip of the jth 
spin. Using Eq. we flnd 

From Eq. (fT^ we obtain the detuning 

A,(g,0) = V-^o-^. = ^E |^_^.|3 =^E(— 7)1- (15) 
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We note that if the long-range dipole-dipole interaction is not taken into consideration, 
then the error to the probabihty amphtude generated by each pulse due to the next-nearest 
neighbor interaction is of the order of (1/2)^=1/8. Our approach allows to compensate the 
unwanted effect of the long-range interaction by optimal choice of the parameters of the 
pulses. 

V. ERRORS AND z COMPONENT OF THE MAGNETIZATION 

The protocol consisting of rectangular pulses with the parameters defined by Eqs. ()12|1 - 
fll5|) allows one to implement entanglement in the ensemble of noninteracting spin chains 
with minimum possible error P. It is possible to relate the error P with the z component 
of the magnetization of the system. Since is caused by the error states, then one 
can assume that A^^ ~ P. We define the dimensionless z component of the magnetization 
M ~ in the following way. Let the maximum value of \M\ be 1. Then for the state 
|00 . . . 00) we have M = 1 and for the state |11 . . . 11) we have M = —1. The measured 
z-component of the magnetization of the system shown in Fig. 1 corresponding to M is 
= \MnRL. 

In this Section we will estimate and compare contribution of different kinds of errors to 
the total error P and relate these errors with M. Since we can estimate P analytically 
for a large number of qubits, the relation between M and P allows us to estimate M 
for L ^ 1. Different kinds of errors can generate positive and negative contributions to 
M so that they can cancel each other. For example, the state Cp(|000111) + |111000)) 
contributes to P but does not contribute to M. In spite of this we assume that for a 
definite range of parameters the relation M = gP holds [see Eq. ()24|1 below], where the 
coefficient g is the "geometrical factor" and P is the probability of error. For example, 
if a single spin in the state (l/-\/2)|0000) of superposition ((T)) is flipped down with the 
probability P then g = —1/L = —1/4. If two spins are flipped down with the probability 
P then g = —2/L = —1/2. If a single spin in the state (l/-\/2)|llll) is flipped up then 
g = 1/L = 1/4 and so on. If, for example, different spins in the state (l/v^)|0000) have 
different probabilities Pj to be flipped then M = {—l/L)J2^^QPj [see Eqs. fl46|l and (j5T|l 
below] . 

The probability of error P depends on parameters and number of qubits L. This de- 
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pendence is different for different kinds of errors. For example, for nonresonant transitions 
P L [see Eq. ()2Up below] . By experimental measurement of M for different parameters one 
can define the most important mechanism responsible for the errors. Using this information 
one can optimize the design and parameters to decrease the error and to prepare a quantum 
computer for implementation of complex quantum algorithms. 



A. Decoherence 



The influence of environment can be characterized by the temperature-dependent spin- 
lattice relaxation time Ti and spin-spin relaxation time T2, where T2 < Ti. The relaxation 
time T2 defines the maximum number of pulses and consequently the maximum number Lmax 
of qubits in each spin chain. Increasing Lmax is desirable for increasing the number of qubits 
in the quantum register and for better measurement. (Because the total magnetization 
of the ensemble of the spin chains is proportional to Lmax-) We note that the number of 
the spin chains -Rmax theoretically is not limited, so that the size of the whole system and 
the total number of spins can be increased by increasing R. If L ^ -/^max the influence of 
environment is small and one can formulate quantum dynamics in terms of wave function 
instead of using density matrix. This allows one to consider quantum logic operations with 



many qubits 



25 



28| and to analytically estimate the influence of other sources of error 
which, as shown below, can cause a much more profound destructive effect on quantum 
computation. 

We define Lmax from the condition 

rH+ E ^3<T2. (16) 
i=i 

where Tj = n/Qj. The pulse performing the Hadamard transform can be made very short 
by increasing the amplitude Bjj of this pulse, so that we will neglect contribution of th- 
Using Eqs. (fTij) - (fTBj) we obtain 

" < 1;^!^^ (17) 

The left-hand side for L^ax ^ 1 can be approximated as 

-1 

^'"^'^ ~ ^ - 0.3399, (18) 





C(3) 
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where ({3) ~ 1.202057, and ({x) is the Riemann zeta function. 

We now estimate Lmax for endohedral fullerenes. The very sharp electron spin resonance 
spectra from the endohedrals ^^N@C6o and ^^P@C6o indicates very long longitudinal relax- 
ation time Ti ~ 10 ms at the temperature T = 10 K, and Ti ~ 0.9 ms at T = 300 K. 
The transverse relaxation time is T2 ~ 20 yus at T = 50 K and T2 ~ 13 /is at T = 300 
K [2^. No nuclear relaxation times have yet been recorded but they are expected to be 
several orders of magnitude longer than the electronic relaxation times. We will make our 
estimation of Lmax for the electron spins. The value of the coupling constant |Je| for two 
electron spins is (/ib = 9.274 x 10^^^ erg/G) |Je|/(27r) ^ 52 MHz. The minimum distance 
between two endohedrals is close to the diameter of the Ceo cage (1 nm) and can be as small 
as 1.1 nm l^. Taking K = 1, T2 = 20 /is, and A = 2.2, the right-hand side of Eq. ()17|1 



becomes equal to 113. From Eq. (fTH|) this value corresponds to Lmax = 136. The value of 
Lmax decreases when the distance A between neighboring endohedrals increases. 

Due to the interaction of the spin system with the environment M becomes positive. The 
reason is that the energy of the state |1) is larger than the energy of the state |0), so that the 
influence of the environment causes the transitions |1) — > |0) while the transitions |0) |1) 
are suppressed. 

B. Nonresonant transitions 

Since wavelength of the radio-frequency pulses is much larger than the size of the quantum 
register, the pulses affect all spins. If the pulse frequency is close to the frequency Uj of jth 



spin then the probability of flipping kth spin, k ^ j is of order of 24 1 

\ 2 



^ \2\j-k\6ujj ■ 

The probabilities of the nonresonant transitions are small provided that the ratio Qj/{2\6uj\) 
is small. Using Eqs. (fT^ and (fl3j) we obtain 

\J\ 

'^'^Aij-kr' ~ V^K^A^\6u;\- ^ ^ 

The probabilities of unwanted quantum states created in the result of the nonresonant tran- 
sitions are of the order of a^, and the probability error Pm- caused by the nonselective excita- 
tions is proportional to a^. The error Pm- grows linearly with the number of pulses L Q. A 
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FIG. 2: The probability error i-*nr obtained using numerical solution 25] and estimate using 
Eqs. (|2U|) and 1)21(1 as a function of the number of qubits L for two values of a. 

typical behavior of Pnr{L) is shown in Fig.|21for a small number of qubits L and for two values 
of a: a = 0.02 and a = 0.09. For example, for two electron spins with |J|/(27r) = 52 MHz 
separated by the distance 2.2 nm and for K=l these values of a correspond, respectively, to 
|5ti;|/(27r) = 141 MHz and |5a;|/(27r) = 31.4 MHz. These values of 6uj correspond, respec- 
tively, to the following magnetic field gradients {6uj = 'je^B^, 7e/(2vr) ^ 28.025 GHz/T): 
2.3 X 10® T/m and 5.1 x 10^ T/m, which can be realized experimentally 30l ISll l32l l33l|. 
We approximate Pnr as 

P„,(L) = -P° +P,\L, L>2. (20) 

The values of P°j. and P^j. found numerically are shown in Fig. IHl The numerical simulations 
are performed by diagonalization of the full Hamiltonian Hn given by Eq. (jSJ for each nth 
pulse in the rotating frame where if„ is time- independent. The obtained eigenstates were 
used for simulation of the quantum dynamics 
presented in Fig. El gives us the expressions 



The best fit obtained from the data 



P„° = 0.8236a^■^^^ P„\ = 0.8615a^-^^^ a < 1 



(21) 



In Fig. 121 we use these parameters to approximate the function P^^{L). One can see that 
there is a good correspondence between the numerical results and our estimate (j20|l . Our 
estimate ()2())1 is especially useful when L is large, L ~ 10^ — 10'^, when no exact solution is 
available. 
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FIG. 3: and P^, in Eq. ^ obtained using numerical solution and the best fits with the 

parameters defined by Eq. (PT|) . 
0.000 1 
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L L 

FIG. 4: The dimensionless z component Mnr of the magnetization generated in result of the 
nonresonant transitions as a function of the number of qubits L for two values of a. The estimates 
are calculated using Eq. I\22^ with the parameters defined by Eq. ()23|). 

Since the z component of the magnetization Mm- is proportional to Pnr it is reasonable to 
present the dependence of Mm- on L in the form 

= M°, - M„\L. (22) 

A typical behavior of M^r as a function of L is shown in Fig. |3 for two values of a. From 
this figure one can see that M^j. is negative. 

From Fig. HI one can numerically calculate M°j, and M^^. We calculated Mq and Mi for 
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FIG. 5: M^j, and M^j. as a function of a. The best fit gives the parameters defined by Eq. ((^HJ. 



different values of a (see Fig. and obtained 



M° = 1.341a2-°^^ Ml = 0.60786a 



1.9795 



a < 1. 



(23) 



From Fig.Elone can see that Eq. ()22p with parameters ()23p gives us a good approximation of 
Mnr. It is important that the number of qubits in Eqs. (j20|) and (j22|) is an exphcit parameter, 
so that one can calculate the probability errors and z component of the magnetization due 
to the nonresonant transitions for an arbitrary number of qubits L. Combining Eqs. ()20|) 
and we obtain the relation between M„r and Pm 



-^nr fi'nr ~^ fl'nr-^nr; 



(24) 



where 



9° 

i»nr 



pO 

nr p\ nr' 



^nr 



Ml 



nr 



C. Interaction between different chains 



Here we will discuss the possibility to decrease the influence of the chains on each other 
by optimal choice of the parameters of the pulses. We will show that if the long-range 
interaction between the chains is not taken into consideration the influence of the chains on 
each other causes the error of the order of Liajd)^ . Consider the jth spin of the rth spin 
chain in the field of the spins of the (r — l)th, rth and (r + l)th chains. (See Fig. IHl) The 
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resonant frequency is [compare with Eq. ()13p ] 



2J 



(25) 



where D is the dimensionless distance between neighboring chains measured in nanometers 
and x = 0'/d = A/D<^l. The Rabi frequency of the pulse is [see Eqs. (fT^ and (fT3|) ] 
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FIG. 6: A scheme for calculation of the influence of the /th spin of the {k — l)th and (k + l)th spin 
chains on the jth spin of the fcth spin chain. 

We now can estimate the error caused by the influence of different chains on each other 
if we neglect the interaction between them. Consider the resonant transition for the excited 
state of the superposition using Eq. (fTUj). If we disregard the influence of the neighboring 
chains on each other, then instead of the resonant transition with Aj(g',p') = we have the 
nonzero contribution of the third term in Eq. (j^^ 



|A,(g',p')l 



2M1 



L-l 



Since the value of the sine in the first Eq. ()1U|) is of the order of unity and the value of the 
cosine is close to zero, the value of the coefficient Cpitj + r^) is of the order of (instead of 
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zero in the ideal case) 



\Cp\t,+T,)\ 



A,(g',p') 













L-l 



E 



slip') 



t [l+X^(j-/)2]3/2 



(27) 



The error Pint caused by the influence of neighboring chains on each other is of the order of 
-Pint = \Cp'(tj + From Eq. ^T7\ one can see that Pint increases with the Rabi frequency 
Qj decreasing. 

In order to estimate the error in the probability amplitude given by Eq. ()27|) consider a 
typical example. Let the distance between the chains be equal to the size of a single chain, 



d = (L - l)a [D = (L - 1)A], L > 1. Then x = " !)• Let us estimate for 
j = L — 1. The first sum in Eq. (jTfjl is 

-3/2 



1 L-2 

E 



1 



1 



2;^^ [l + x2(L-l-/)2]3/2 
For all terms of the sum we have 

1 
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E 

1=0 



1 



I 



L- 1 



23/2 



< 



1+1 



L-l 



-3/2 



< 1, 



SO that 



23/2 

For the second sum we have 

1 



L-l ^2 
< E 

1=0 



1+1 



L- 1 



-3/2 



<L-1. 



L-2 

E 

1=0 



1 1 

l + :^ + :^ + --- + 



(L - 1 - /)3 ' 23 ' 33 
Finally, assuming K = 1 we obtain 



(L-iy 



< C(3) ^ 1.202. 



Pint>0.51(L-l)x^ 



(28) 



This is the error introduced by only one pulse. The errors generated by different pulses of 
the protocol can accumulate. 

From Eq. ()28|1 one can see that a small parameter which characterizes the dipole-dipole 
interaction between the chains is Lx^ rather than x^- The influence of neighboring chains 
on each other can be minimized by introducing the corrections to the frequency [the third 
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term in the right-hand side of Eq. ()25|) ] and to the Rabi frequency [the second term in the 
right-hand side of Eq. (f^Bj) ]. These corrections minimize the errors only for intermediate 
chains and do not minimize the error for the edge chains with r = and r = R — 1. One 
can use our corrections if the number of chains R is large, so that one can neglect the edge 
chains, or when the chains are placed relatively close to each other, when these corrections 
are relatively large. 

If one minimizes the errors caused by the nearest-neighboring chains, the error in the 
probability amplitude caused by the influence of the next-nearest neighboring chains is of 
the order of \/Pmt/8. Using our approach one can minimize the errors caused by the next- 
nearest neighboring chains and disregard the chains with r = 0, l,i? — 2,i? — 1. This makes 
the error caused by the next-next-nearest neighboring chains to be of the order of v^^nt/27 
and so on. 



D. Qubit displacements 

One of the most serious problems that prevents building a solid-state quantum computer 
is manufacturing the spin system like that shown in Fig. 1. Atoms with nonzero spin, such 
as ^^P, can be placed on the surface of a magnetically neutral substance, such as ^^Si, using, 
for example, scanning tunneling microscopy technique The placement of the qubits 

can be not perfect, so that these qubits form distorted spin chains. If one deals with a 
single chain, one can measure the locations of the qubits and to chose the suitable pulse 
parameters to compensate the deviations of the qubits from their prescribed positions. On 
the other hand, if one implements a quantum algorithm on an ensemble of spin chains, the 
deviation in the location of a qubit from the perfect position (called below displacement) in 
a chain makes this chain different from other chains, and this error cannot be completely 
compensated by a proper choice of the parameters of the pulses. Here we will investigate 
this kind of error. 

Since a qubit in a solid state is usually incorporated into the crystal lattice, the minimum 
possible qubit displacement is equal to the lattice constant. If the displacement happens 
in the direction of the magnetic field gradient (along the x axis in Figs. 1 andlH)), then 
even a small displacement causes a relatively large change in the Larmor frequency of this 
qubit because the magnetic field gradient is supposed to be large. Consequently, we believe 
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that this kind of error causes the most profound destructive effect on quantum computation 
in our system. Since the frequency of the displaced kth qubit considerably differs from the 
frequency of the pulse, this qubit will not flip and the other (/c + l)th, (/c + 2)th, . . . , (L — l)th 
qubits will not flip also. For example, if k = 3 and L = 6 the excited state evolves as 

IO5O4O3O2Q1I0) ^ lOOOOll) ^ lOOOOUl) ^ lOOOlll) ^ loooiu) ^ loooiu), 

where the qubit to be flipped by the corresponding pulse is underlined. One can see that 
the z component M of the magnetization due to the error caused by displacement of qubits 
is positive. 

We now calculate the error due to the qubit displacement (s). It is convenient to define 
the dimensionless displacement Vk of the kth qubit as 

= 4^ = (29) 

duj a 

where dx^ is the dimensional displacement, dui^ is the change of the Larmor frequency of 
the kth qubit caused by this displacement, the sign "+" must be used if one considers 
nuclear spins and the sign "-" must be used for electron spins. For example, the value 
Vk = 1/15 corresponds to the displacement by one lattice site if the number of atoms 
between neighboring qubits is equal to 14 and by two lattice sites if the number of atoms 
between neighboring qubits is equal to 29. 

If the probability V of a. displacement is relatively large, V > 1/L then the number 
of "perfect" spin chains, where all spins are not displaced, is relatively small, so that it 
is important first to study the errors Pd and the magnetization Md caused by the qubit 
displacements in a single spin chain. 

The displaced kth qubit affects all other qubits in the chain. The transition frequency 
and detuning of the jth qubit (j 7^ k) change by the value [see Eq. (fT^ ] 
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< nj. (30) 



It is convenient to characterize the influence of the displaced qubit on all other qubits by a 
small dimensionless parameter 



_ SvkV^K^ - 1 
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< 1. (31) 
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Then in Eq. fjlO|) we have 



When the jth qubit is flipped in the excited state we have A = Aj^ instead of A = and 
the error is 



For example, for = 1/20 and K = 1 we have 



Pk±i,k ~ 0.006, Pk±2,k ~ 2.3 X 10- 



The probability of the transition from the ground state caused by the deviation in the 
detuning is 



For K = 1 the error caused by the ground state Pjj^, ~ 0.25/3|j;. is of the same order as the 
error Pj^ caused by the excited state. Below we will neglect the errors Pjk and Pj^ caused 
by the influence of the displaced qubit on all other qubits as being small compared to the 
other errors. 

We now estimate the probability Pkk caused by the kth pulse on the kth displaced qubit. 
Since a ~ l/|<^t^| it is convenient to measure the dimensionless frequency displacement in 
units of I /a, so that the change in the transition frequency of the kth qubit caused by 
its displacement is proportional to f^/a. We flrst analyze the action of the A;th pulse on 
the excited state. Instead of the resonant transition with the detuning A^ = we have 
the transition with the detuning A^ = ±Vk6uj, where the sign "+" corresponds to the 
displacement in the positive x direction and the sign "-" correspond to the displacement in 
the opposite direction. 

If /c 7^ then the amplitude of the excited state is 



\Cq>\ = —j= sni 

V2 ^1 + {vkja^f 



2V Wj 



where 
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The error generated by the excited state is 



kk 



1 



sm 



Vk\' 
Oik) 



(33) 



+ (vk/ak) 

One can see that the ratio Vk/a characterizes the error caused by the displacement Vk- if 
Vk/c( — > then 0, otherwise 1/2. 

The error caused by the action of the kth pulse, /c 7^ 0, on the ground state of the 
superposition is 



1 



1 + (V^K^ - 1 T ^ 



sm 



(34) 



For a large magnetic field gradient the error P^ is small 

1 w 



> 1. 



(35) 



{vk/akf «fc 

For A; = the probability of the excited state after implementation of the Hadamard 
transform on the displaced qubit is 



1 + {vo/aQY 



■ sm 



(36) 



where = VLij/\5uj\, Qh is the Rabi frequency of the pulse implementing the Hadamard 
transform. Here and in the sequel we take 0:0 = 0:1= o. The probability of error generated 
by this pulse is 



+ 



\Ci\ 



where |CoP = 1 — jC'ip. Assuming that the other pulses of the protocol do not generate 
error, we obtain that the error due to the displaced zeroth qubit is 



P 



00 



P" 



(1-r/) 



+ 



7] 



1 - 2r/. 



(37) 



Since l/o ~ l^c^l decreasing \5uj\ decrease the error P^. On the other hand, due to Eq. ((201) 
and ()21|). the error P^^ caused by the nonselective excitations (nonresonant transitions) 
increases with \5uj\ decreasing. The total probability error for 7^ is 



P = Pr^r + Pa + P'A 



(38) 
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FIG. 7: The probability error P obtained using numerical solution as a function of 1/a for different 
values of Vk- The displaced qubit is located at the center of the chain, L = 9. 

where the qubit is assumed to be displaced in the negative x direction. In Fig. [7| we plot 
the probability error P, defined as 

1 



-\CoiT)[ 



+ 



\C2^-iiT)\' 



(39) 



which was obtained using exact numerical solution for A; 7^ 0. In Eq. (j39j) T is the total 
time of implementation of the entanglement protocol. When l/a is small the probability 
error P is large due to the nonresonant excitations. When 1/a is large P is large because 
the displaced qubit does not flip. From the results presented in Fig. [7| one can see that 
if the displacement is relatively large {vk = 1/10), then the error is always large and the 
entanglement protocol cannot be optimized for any parameters of the model. 

As follows from Fig. [7[ the magnitude of P becomes close to 1/2 and relatively indepen- 
dent of L and 1/a when 1/a is large, i.e. when (f^/a)^ ^ 1. For these parameters we have 
mostly two states in the superposition: the ground state and the partially excited state. 
For example, if the kth qubit is displaced, k = 3 and L = 6, then instead of the desired 
entangled state 



V2 



(looooo) + liiiiii)) 



(40) 
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we have the state 

^ (lOOOOO) + lOOOlll)). 
v2 

The z component of the magnetization due to a displaced qubit can be estimated 
using the probabihty error. If 7^ the probabihty is mostly independent of the position 
k of the displaced qubit in the chain. Unlike P, Md is large (and positive) if the displaced 
qubit is located in the beginning of the spin chain and relatively small if the displaced qubit 
is located in the end of the chain. For example, if k is the number of the displaced qubit 
and k = \ then the excited state is (L = 6) |000001) instead of |111111) and the entangled 
state is 

^ (lOOOOO) + lOOOOOl)) (41) 
v2 

instead of the state (001) • The z component of the magnetization for to the state (PT|) is 
= (5/6). If the displaced qubit is located in the end of the chain, for example, if /c = 6, 
then the value of M for the state 

^ (lOOOOO) + lOlllll)) (42) 
V2 

is Md = (1/6). 

In order to relate M with the probability error consider the two situations when 7^ 
and when = 0. If /c 7^ we note that after implementation of the entanglement protocol on 
the spin chain with a displaced qubit in the register there are mostly four quantum states: 
the ground state with the probability 1/2 — P^(±), the error state with the probability 
Pd(=t) created from the ground state, the error state with the probability Pd created from 
the excited state, and the fully excited state |11 . . . 11) with the probability 1/2 — Pd. Next, 
we assume that the position of the displaced qubit in the chain is random. By averaging 
over many random realizations we obtain that the two error states do not contribute to M. 
For example, the z component of the magnetization of the partially excited state in Eq. ()4ip 
is M = I ■ I while in Eq. ()42|1 M = — | ■ |, so that the average of these two contributions is 
zero. The contribution to the M due to the fully excited state is — Q — Pd^. By adding all 
these contributions we obtain 

Md = ^ - P^(±) - - Pd) = Pd - P'^{±). (43) 

If A; = there are mostly two states in the register: the ground state with the probability 
1 — r/ and the fully excited state with the probability 77, where 77 is given by Eq. (jSHl)- The 
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FIG. 8: The z component of the magnetization M as a function of 1/a for R = 100 spin chains 
with L = 7 qubits in each chain. The numerical results are averaged over 50 realizations of different 
7x 100 qubit ensembles with randomly chosen displaced qubits, ^ = 1/(5L). The displacements are 
in the random directions along the x axis. The estimate is calculated using Eq. (|47l) . Vk = 1/20, 
K = 1. 

magnetization due to the displaced zeroth qubit is 

= 1 - 2r/ = PH. (44) 

The total magnetization is 

M = M^, + Md, for ky^O, 

M = M^, + M'H, for A; = 0, (45) 

where M^^ is given by Eqs. ((221) and (j^ . 

Assume that we have an ensemble of R spin chains and the probability of a qubit to be 
displaced by one lattice site is ^. If, for example, ^ = 1//^, then on average one qubit in each 
chain is displaced, if ^ = l/{nL) then on average one qubit in n chains is displaced. The 
total number of displaced qubits with /c 7^ is on average ^R{L — 1) and the total number 
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of displaced qubits with /c = is on average ^R, so that 

pu^lnkH+pu+)], (46) 

where Pd,fc and P^^ depend on k through the dependence of ak on k in Eq. (jH21)- We 
calculated numerically 

M = Mnr + Md (47) 

as a quantum-mechanical average for 100 noninteracting spin chains with randomly chosen 
displaced qubits and random displacement directions (in positive and negative directions 
along the x axis). In Fig. |Hlwe plot M as a function of a obtained using the exact numerical 
solution and our estimate given by Eq. (jTfj) for ^ = 1/ (5L) . As follows from the figure our 
estimate (jlTj) is a good analytical approximation of M. In particular, Eq. (P7j) can be used 
for estimation of M when the number of qubits L in each chain is large. 
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E. Fluctuations of permanent magnetic field 



The error Pose caused by fluctuations of permanent magnetic field duo due, for example, 
to unwanted oscillations of the DC current in wires needed for creating the magnetic field 
gradient can be estimated using Eqs. and (jH^ . Instead of the dimensionless deviation 
Vk we introduce the average dimensionless deviation v as 



5uj 



where 



is the average deviation of the transition frequency of a qubit from the optimal 
value caused by a fluctuation of the current in the wires creating the magnetic field gradient. 
We assume that this deviation is small, {v/aY ^ 1- 

The average error due to unwanted transitions from the excited and ground states for 
/c 7^ is [see Eqs. (jH^ and 



Posc{k) — Pogc(A;) + Pogj,(A;), 



(48) 



where 
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Here -Posc(^) ^oscW probability errors created from, respectively, the ground 

and excited states by action of the fcth pulse of the protocol. For k = (Hadamard gate) 
from Eq. the probability of the excited state is 



V 



1 + 



so that the average error is 

Poscik = 0) 



2r] 



1 



V \2 

aoJ 



The average total error due to the nonresonant transitions and oscillations of the magnetic 
field is 

L-l 

P = P,r + Posc(0) + E^osc(fc) 



(49) 



k=l 



Since P„ 



Q^opt 



a 



and Posc{k) ~ there is an optimal value of 
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-0.82 + 0.86L) 



1/4 
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TT 
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1 7i\4K^ - 

2 ^ 



L 



e{3) 



+ 0.6 



1/4 



(50) 



where P is minimal. For example, for v = 10~^, L = 9, and = 1 we have aopt ~ 9.08 x 10^^. 

In Fig. El we plot the error P caused by action of both nonresonant transitions and 
unwanted oscillations of the permanent magnetic field. Numerical results are obtained by 
diagonalization of the full Hamiltonian matrix in the rotating frame. The random values 
of V for the numerical results have the Gaussian distribution with the zero average and 
dispersion equal to v. In our model a magnetic field fluctuation is constant (and random) 
during each pulse. Each point of the numerical results is the average over 30 realizations 
with the random values of v. One can see from the figure that for a < Oopt the error is 
mostly defined by the fluctuations of the permanent magnetic field and P decreases with 
increasing a as P ~ l/a^. For a > Oopt the error is due to the nonresonant excitations and 
P ~ a^. From Fig. El one can see that our estimate correctly describes the probability 
error in the presence of unwanted oscillations of the permanent magnetic field. 

Now we will calculate the z component of the magnetization Mqsc- Consider a typical 
example with L = 4 qubits. If the first pulse of the protocol {k = 1) generates error and other 
pulses do not generate the error, then after implementation of the entanglement algorithm 
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FIG. 9: The error caused by nonresonant transitions and unwanted oscillations of the permanent 
magnetic field as a function of alpha for v = 10~^ and L = 9. The estimate is obtained using 
Eq. (EH). 
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FIG. 10: M as a function of a in the presence of unwanted oscillations of the external permanent 
magnetic field. The parameters are the same as in Fig. El The estimate is calculated using Eq. (|52() . 
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there are four states in the register with the following probabilities: 

10000) , l-PUl), 

11110) , P4,(i), 

11111) , 1-p:M) 

10001) , P^Jl). 
The value of Mqsc for this superposition is 

If all the pulses of the protocol generate the error, then the z component of the magnetization 
is 

L-l 

Mosc = Posc(O) + E 
k=l 

In Fig. ^1 we compare our estimate 

M = + Mo,e (52) 

with the results of numerical simulations for the same values of a as in Fig. 121 One can see 
that our estimate gives us a good approximation of M. 

VI. DISCUSSION 

We considered the implementation of entanglement in a two-dimensional ensemble of spin 
chains. We demonstrated that the entanglement can be created in the system with a long- 
range interaction, such as dipole-dipole interaction, and the error caused by this interaction 
can be significantly reduced by optimization of parameters of the pulses. 

If the entanglement is implemented with no error, then M = 0. We considered different 
mechanisms which can generate errors and make M ^ 0. By experimental measurement of 
M for different parameters one can define the most important mechanism responsible for 
the errors and optimize the design using the obtained information. The output signal can 
be enhanced by the multiple copies of the spin states. The most important dimensionless 
parameter characterizing the model is a which is proportional to the ratio of the Rabi 
frequency to the difference 6uj between the Larmor frequencies of neighboring qubits. We 
now summarize the influence of different kinds of errors on M. 



1 + 



L~2k 
~L 



[PLik) - PUk)] . 



(51) 
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1. Decoherence. Decreasing a by decreasing the Rabi frequency one can increase the 
total time of implementation of the algorithm and increase |M|. The value of M is 
positive. Decreasing a by increasing the gradient does not influence M. 

2. Nonselective excitations (nonresonant transitions). Decreasing a by decreasing the 
Rabi frequency or increasing one can decrease |M|. The value of M is negative. 

3. Influence of different chains on each other can be estimated in the following way. (a) 
One implements the protocol using the frequency v and the Rabi frequency f2 given by 
Eqs. fll3|) - (fT3j) . (b) One takes into consideration the influence of neighboring chains on 
each other and modifies the frequencies and the Rabi frequencies using Eqs. and 

If |M| decreases then the error is mostly caused by the dipole-dipole interactions 
between the spins of different chains. 

4. The error is caused by displacements of the qubits if one observes the following prop- 
erties: (a) if the gradient buj is relatively large [(ffc/a)^ 1 in Eq. (jHSj)] is positive 
and independent of a; (b) the same effect is observed if the distance a between the 
qubits is relatively small and the relative displacement of the fcth qubit \v}\ = \dxk\/a, 
where dxk is the displacement of the kth qubit along the chain, is relatively large, 
{vkl > 1/10; (c) if the gradient 6uj is relatively small (|ffc|/a > 1) increasing a we de- 
crease \M\ and M > 0. The latter effect is opposite to the influence of a in nonresonant 
excitations where increasing a we increase \M\ and M < 0. 

5. Unwanted fluctuations of permanent magnetic field. Increasing a decrease M which 
is positive and M ~ 

We did not consider all possible mechanisms such as, for example, the influence of mag- 
netic impurities in the substrate 3]. Some processes not considered here in detail, such as 
the decoherence caused by environment, can be investigated by using the density matrix if 
one finds that this kind of decoherence is the most important mechanism responsible for the 
errors. 
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